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Abstract 

We revisit electromagnetic field propagation through tight-binding arrays of coupled photonic 
waveguides, with properties independent of the propagation distance, and recast it as a symmetry 
problem. We focus our analysis on photonic lattices with underlying symmetries given by three 
well-known groups, SU( 2), S77(l,l) and Heisenberg-Weyl, to show that disperssion relations, 
normal states and impulse functions can be constructed following a Gilmore-Perelomov coherent 
state approach. Furthermore, this symmetry based approach can be followed for each an every 
lattice with an underlying symmetry given by a dynamical group. 
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I. INTRODUCTION 


Photonic lattices, that is, tight-binding arrays of photonic waveguides, have proved a 
valuable tool in the classical simulation of quantum physics [1. 2). Vice versa, the insight 
gained from these quantum-classical analogies has been useful as a tool for photonic lat¬ 
tice design [3|. This comes from the fact that the differential equation set describing the 
propagation of held amplitudes through an array of nearest neighbor coupled waveguides 

SHE], 

- = uf(l, z)Si(z) + A [g(l, z)S i - 1 (z) + g(l + 1, z)S t+1 (z )], (1) 

with / = 0,1, 2,... and the restriction g(0, z) = 0, can be cast in a vector differential equation 
form similar to Schrodinger equation, 


~iX\S(z)) = H(z)\£(z)). 


( 2 ) 


Here, the held at the Z-th waveguide is £i(z), the effective waveguide refractive index is 
uif(l,z), and the effective nearest neighbor coupling A g(l,z) [3]. We have borrowed Dirac 
notation from quantum mechanics, such that vectors are written as kets, 


\S(z)) = 


( So (z)^ 
Si(z) 

£ 2 (z) 


V 


7 


^2 £l (z)\j), 


( 3 ) 


( 4 ) 


1=0 


that collect all the information from the held amplitudes, and the real-symmetric tridiagonal 
matrices that collect all the properties of the tight binding array are written as operators, 


/ 


H{z) = 


0,z) Xg(0,z) 0 

A^(0,*) uf(l,z) Xg(l,z) 

0 Xg(l,z) uf(2,z) Xg(2,z) 


\ 


V 


7 


= uiA.q^z') + A A+^z) + A_( y z) 


( 5 ) 


( 6 ) 


The operator A$(z) = f(n,z ) is a diagonal matrix, the operators A + (z) = g(n,z)V t and 
A_(z) = Vg(n,z) are lower- and upper-diagonal matrices in terms of the step, n\j) = j\j), 
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up-step, V ] \j) = |j + 1), and down-step, V\j) = |j — 1) operators. Note that the functions 
related to the effective index and coupling parameters, f(l, z ) and g(l, z), have been rewritten 
in terms of the step operator, with action f(n,z)\l) = f(l,z)\l) and g(h,z)\l) = g(l,z)\l), 
such that it is straightforward to recover the original differential set. At this point, it is 
clear that the held amplitude distribution at a point z in the photonic lattices is given by 
the solution to the Schrodinger-like equation in (J2|, 

\£{z)) = e i % 6{ - x)dx \£(ti)), (7) 


where the information from the initial held amplitudes distribution impinging the lattice is 
collected in the vector |£(0)). It is straightforward to connect this result with the impulse 
function, T myn (uj , A, z ), that provides the held amplitude at the n-th waveguide given that 
the initial held impinged just the m-th waveguide, 


= (n\e if °* ix)dx \m). (8) 

Note that it is typically hard to calculate the propagator e* -h H ( x ) dx but underlying symme¬ 
tries may provide help mm- 

Our goal here is to bring forward the fact that underlying symmetries may simplify greatly 
the calculations of dispersion relations, normal modes and impulse functions for photonic 
lattices. We will focus in the simplest case of lattices which parameters do not depend 
on the propagation distance. In this case, if the operators A 0 and A± form a closed group 


under commutation 


A, B = AB — BA, it is possible to use so-called disentangling formulas 


TTfl [To] to calculate the corresponding impulse function, 


l min (u,X,z) = {n\e tzH \m), 

= (n|e i2 H 0+A ( i++i -)]|m), 

_ / n \ e a+(z)A+ e ao(z)Ao e a-(z)A- | m ) 


(9) 

( 10 ) 

(ii) 


in terms of a particular class of generalized Gilmore-Perelomov coherent states [I61 - H8] . 

D [a 0 (z),a ± (z)] \m) = e a+W^+ e aoW^ e a_(*)A_ | m )_ ( 12 ) 


Thus, these photonic lattices can be seen as optical simulators of generalized Gilmore- 
Perelomov coherent states. In the following sections, we will provide three working ex¬ 
amples of lattices with an underlying symmetry and show that a symmetry based, Gilmore- 
Perelomov generalized coherent state approach can help us providing dispersion relations, 
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normal modes and impulse functions. First we will work with finite lattices with underly¬ 
ing 577(2) symmetry; the so-called J x photonic lattices used to produce coherent transfer 
between input and output waveguides mm are a particular example of this underlying 
symmetry. Then we will move to semi-infinite lattices described by the 517(1,1) group which 
have as particular realization in the literature the so-called Jacobi lattices (2U. Finally, we 
will explore lattices with underlying Heisenberg-Weyl symmetry that includes as a specific 
example the so-called Glauber-Fock lattices [221126] and close with a summary. In all the 
following cases, the closed forms for the respective impulse functions, in terms of the cor¬ 
responding generalized Gilmore-Perelomov coherent states, are new in both the optics and 
quantum mechanics literature, to the best of our knowledge. In order to corroborate them, 
all the analytic results were compared with numeric diagonalization, for normal modes, 
and with results from a normal-modes-decomposition and small-step propagation, for the 
impulse functions. In all cases analytic and numeric results are in good agreement. 


II. LATTICES WITH UNDERLYING SU( 2) SYMMETRY. 


Let us start with a finite array of coupled waveguides described by the matrix, 


H S u( 2 )(u, A,j) =u(n- j) + A V \/n(2j + 1 - h) + yj h(2j + 1 - h)V ] 
= cuJq + A(J_ + </_)_), 


(13) 

(14) 


in terms of the 517(2) group [j + , j_] = 2 j 0 and [J 0 , j±] = ± J± [27] and j a positive integer 
such that 2 j + 1 is the size of the lattice. It is straightforward to construct a Gilmore- 
Perelomov displacement operator mm, 


Dsu( 2){Q) — e 


_ _ 


__ e jfy tan|e|J +e l nsec 2 \ 0 \ Jog-pj tan \ 6 \ J_ 


(15) 

(16) 


that diagonalizes this matrix, 


Dsm(0)Hsu(2)(u,\)Dl u{2) (e) = Vuj 2 + 4A 2 (n - j), 


with 


2A 

tan 2 6 = —. 

u 


(17) 


( 18 ) 
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Thus, we can immediately construct the dispersion relation for this photonic lattice, 


Msu( 2 ) (m, uj, A, j) = Vw 2 + 4A 2 (m - j ), (19) 

and its normal modes, given by what we are going to call 577(2) displaced number states, 
such that the n-th component of the m-th normal mode is given by the expression, 

(n\n S u( 2 )(m,uj,X : j)) = (n\D SU (2){-0)\m), (20) 


\ 


2 j \ I 2j 

m I \ n 


= (- 1 ) 

xK m (n, sin 2 \0\,2j) , 


(cos|d|) 2i " m " n (sin \ 9\) m+n x 


( 21 ) 


where we have used the fact that the lattice parameters are real and positive, A > 0 and u ; > 


0. The notation 


and K n (x,p, N) stand for the binomial coefficient and Krawtchouk 


polynomials (28], which are well known in the discrete optics community [29H32] . in that 
order. Now, it is also possible to compute the impulse function of this array by using 
standard disentangling formulas for 577(2) [T51 - IT5] . and some algebra, in terms of generalized 
Gilmore-Perelomov coherent states, 


l^\u,X,j,z) = {n\e izA ^\m) } (22) 

/ 2 j\ (2iAsin | \Zu 2 + 4A 2 )” t+ " 

\ n ) (^ 2 + 4A 2 )- 7 X 

/ _ 2 , _ Z _\ 2 j—m—n 

x yVoj 2 + 4A 2 cos - y/u 2 + 4A 2 — ioj sin - Voj 2 + 4A 2 J x 

f A\2 _ _ \ 

xKm \ ’ U 2 + 4A 2 sin2 + 4A2 ’ 2j J • 

Note that the impulse function will show a 2tt/\Juj 2 + 4A 2 periodicity due to the squared sine 
argument in the Krawtchouk polynomials. Figure [l] shows the periodicity discussed above 
in the intensity, that is, the squared impulse function, X,j, z)\ 2 , of an initial 

field amplitude impinging just the central waveguide of an 517(2) lattice with parameters 
A = 0.4a; and j = 5. 

In the specific case provided by waveguides with identical effective refractive index, such 
that the inclusion of an overall constant phase leads to oj — 0, the Hamiltonian-like matrix 
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FIG. 1. Intensity at the n-th waveguide, (co, X, j, z)\ 2 , for light initially impinging the 

middle waveguide, rn = 5, of an 51/(2) array with parameters A = 0.4cu and j = 5. 


becomes, 

Hsu(2)(®i A, j) — 2AJ X , J x = - ^J+ + J- j . (24) 

This may be the origin of the J x lattice moniker for such restriction [EH I2D] • In this case it 
is straightforward to take the limit, 

, 1 2A 7T 

lim - arctan — = — (25) 

2 LU 4 K ' 


and use this value to find the normal modes, 


{n|f!(m,0,A,j) s( 7(2)) = K m (n, 1,2A , 


(26) 


and the impulse function, 


^ (2)) (0,A,i,.) = 


2 j \ / 2 j 


m \ n 


N 

xK m (n, sin 2 Az, 2j) . 


(isin Xz) m+n (cos \z) 2j 171 n x 


(27) 


It is straightforward to notice that the lattice shows a 7 t/A periodicity. It is also simple to 
check the properties of Krawtchouk polynomials and realize that the initial field amplitude 
at the m-th waveguide will be transferred to the (2j— m)-th waveguide at half the periodicity, 
Az = 7 t/2, 

l!n U n™(0A,],^)=S n<2j _ m , (28) 
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where the notation S a ,b stands for Kronecker delta that is equal to one when a = b and zero 
for all other cases. This has been used for coherent transfer in the so-called J x lattices mm 
and is an immediate consequence of the impulse function given in terms of trigonometric 
functions and Krawtchouk polynomials. Figure [2] shows this coherent transfer and the 
periodicity discussed above in the intensity, |Tm'n' 2 ^(ay A, j, z) | 2 , of an initial field amplitude 
impinging just the m — 2 waveguide of an SU (2) lattice with parameters u — 0 and j = 5. 

0 

n 5 

10 


FIG. 2. Intensity at the n-th waveguide, X, j, z)\ 2 , for light initially impinging the 

m = 2 waveguide of an SU( 2) array with parameters uj = 0 and j = 5. 



III. LATTICES WITH UNDERLYING 517(1,1) SYMMETRY 


There is another three-element closed group that can be used to describe photonic lattices, 


Hsu ( 1 , 1 ) (cj, A, k) = uj (h + k) + A V\J n(2k — 1 + h) + y/h(2fc — 1 + h) V * 
= aiA' 0 + A (AT + K+). 


(29) 

(30) 


Here the operators form the SU (1,1) group that satisfy the commutation relations, [AT + , AT] = 
—2A” 0 and [A” 0 , A^] = ±K±. Here, the ideal lattice has infinite size but these results will 
describe a real-world finite lattice as long as the propagated electromagnetic held amplitudes 
stay far from the edge of the array. The real positive parameter k > 0 is the Bargmann 
parameter of the SU(1,1) group [2?J. Again, this effective Hamiltonian-like matrix can be 
diagonalized, 


Dsu(i,i){€)H S u(i,i){u}, A, k)L>l u{1A) (0 = v^ 2 - 4A 2 (n + k), (31) 
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with the difference that a strong restriction appears, 


tanh 2£ = —, uj > 2 A. 
u 

Here the displacement operator has the form, 

Dsu(i,i)(0 = e^ + -^ k -, 

= gjfj tanh I?! K+ ginsech 2 |£|kotanh \£\K- 


(32) 


(33) 

(34) 


where we have used the standard disentangling formulas for 517(1,1) [T314T5] . Under the 
restriction, uj > 2A, the dispersion relation is: 

^517(1,1) (m, u,X,k) — Vuj 2 — 4A 2 (m + k ), (35) 

with the n-th component of the 777 -th normal states given by the 517( 1,1) displaced number 
states, 


(n\^su(i,i)(m,uj,X,k)) = (n\D su{hl) (-^)\m), 


(-!) n T(m + 2k)T(n + 2k) ,^m+n ,^-2fc-m-n 


mini 


r (2 k) 

X 2 Fi (—m, —77, 2k, — cscll 2 |^|) , 


(sinh 1^1) (cosh |^|)' 


(36) 

x 

(37) 


where the notation r(x) and 2 Pi(a, b; c; x) stand for the Gamma and Gauss hypergeometric 
functions, in that order, and we have accounted for the fact that the paremeters uj and A 
are real and uj > 2A. It is quite important to note that, here, it is only possible to calculate 
the normal states for the restriction uj > 2A with the Gilmore-Perelomov coherent approach 
due to the characteristics of the 517(1,1) group [33]. 

Again, it is possible to work out an impulse function in terms of what we will call 517(1,1) 
generalized Gilmore-Perelomov coherent states, 


l^\uj,X,k,z) = (n\e izksu ^\m), 


(38) 


Y{2k) 


(4V - *’). /r(m + 2k)T(n + 2k) £^^y" + " x 


mini 


x 4A 2 — uj 2 cosh -\/4A 2 — w 2 — iuj sinh - \/4A 2 — w 2 ^ 

'uj 2 - 4A 2 


—2 k—m—n 


X 


X 2 -Pi —77, —777 ; 2 k] 


4A 2 


csch 2 - \/4A 2 — w 2 ) . 


(39) 


Here it is of great importance to note that this impulse function works for any given case of 
uj and A as the restriction uj > 2A appears only for the dispersion relation and normal states. 















Thus, we have two types of propagation regimes. The first one is given by the relation oj > 2A 
that transforms the hyperbolic functions into trigonometric functions. Thus, the lattice will 
be 27r/ \]uj 2 — 2A 2 periodical. Figure [ 3 ] shows this periodicity in the squared impulse function, 
A ,j, z)\ 2 , of an initial field amplitude impinging just the m = 15 waveguide of an 
577(1,1) lattice with parameters A = 0.4a; and k — 1/4. 

0 

n 

200 


FIG. 3. Intensity at the n-th waveguide, {oj-, A, k, z)\ 2 , for light initially impinging the 

m = 15 waveguide of a 577(1,1) array with parameters A = 0.4a; and k = 1/4. 

0 

n 

200 

0 — — 

4 2 

Az 

FIG. 4. Intensity at the n-th waveguide, |I^^ 1,:L ^(a;, A, k, z)\ 2 , for light initially impinging the 
m = 15 waveguide of a 517(1,1) array with parameters u = 2A and k = 1/4. 

The other will be defined by oj < 2A, where the lattice is not periodical and light impinging 
the m-th lattice will disperse with propagation due to the hyperbolic functions. In the 
special cases, oj = 2A and oj — 0, the impulse function can be calculated via the limits of 
the hyperbolic functions. The impulse function for a; = 2A, shown in Fig. [4] for Bargmann 
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parameter k = 1/4, is given by the expression, 


i))(2 A ,A,M) = J T ( m + 2k ) T ( n + W (~l) fc+m+ "(A^(A, + x 


m\n\ 


r( 2 fc) 


x 2 -fi ( ——m; 2k; 


X 2 z 2 


(40) 


and for a; = 0, shown in Fig. [5] for k = 1/4, by the expression, 

n-(su(hi))m 4 u _ . /r(m + 2fc)r(n + 2fc) (i sinh Az) m+n (cosh Az) 

(u,A,/c,z)-y ^ mM 


—2 k—m—n 


X 


x 2 -Pi (— n, —m; 2k; —csch 2 Az) . (41) 

Note that a particular realization of lattices with this underlying symmetry described by 
the parameters cj = 1 + /3 2 ,A = /3 with the parameter /3 such that to > 2A for the Bargmann 
parameters k = 1/2 and k — 1, have been used to report pairs of isospectral arrays of 
photonic lattices as an analogy to supersymmetric quantum mechanics. In these cases the 
impulse function reduces to Jacobi polynomials and gave rise to the moniker of Jacobi 
photonic lattices [Tl] . 



FIG. 5. Intensity at the n-th waveguide, A, k, z)\ 2 , for light initially impinging the 

m = 15 waveguide of a 577(1,1) array with parameters to = 0 and k = 1/4. 


IV. LATTICES WITH UNDERLYING HEISENBERG-WEYL SYMMETRY 

The Heisenber-Weyl lattice has an effective coupling that increases as the square root of 
the position of the waveguide, 

H hw (uj, A, z ) = ion + A [y + \fhV^ 

= ion + A(a + a*), 
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(42) 

(43) 











where we have used the creation (annihilation) operators a' (a) that satisfy the commutation 
relations [cd,d] = —1, [n, cd] = a\ and [n,a] = —a. This matrix can be diagonalized, 

A 2 A 

Dhw{oi)HhwD^ hw (oi ) = un -, a = —. (44) 

C 0 C0 

by the Glauber displacement operator [33, 35], 

Dhw{<x) = (45) 


Thus, the dispersion relation, 

A 2 

A) = um -, (46) 

UJ 

and the normal modes are given by the displaced Fock states (3S1 157] . 


(n\tt H w{m,uj, A)) 


{n\D H w(-Oi)\m), 


(-ir.-n* 

\/m!n! 


2 I^I 2 I -- 


U —m, n — m + 1, 


(47) 

(48) 


where the U(a,b,x ) stands for Tricomi confluent hypergeometric function |38j . Typically, 
displaced number states are given in terms of generalized Laguerre polynomilas, but this 
requires to split the result into two cases, one for n > m and another for n < m. We 
favor Tricomi confluent hypergemetric functions because it is straightforward to use them 
in numeric simulations, rather than defining the cases required with generalized Laguerre 
polynomials. 



FIG. 6. Intensity at the n-th waveguide, |A, z)| 2 , for light initially impinging the m = 24 
waveguide of a Heisenberg-Weyl array with parameter A = 0.4 co. 
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The impulse function of this array is given by the following expression, 


Zg^(w,A,z) = (n|e“*«“'|m), 

(- 1 )’” 

= , < 
v mini 


iz I muj — 


jA|2( 1 _ e -^) 


- (e iwz - 1) 

LU V ' 


X 


TT , ^ A |2 . 2 LUZ 

xU ( —m, n — m + 1,4| — | snr — 

LU 2 


(49) 


(50) 


where we have used the properties of Glauber displacement operator [31]. Thus, light im¬ 
pinging just at just the m-th waveguide will produce a field amplitude distribution at the 
waveguides that is equivalent to the amplitude distribution of the quantum optics displaced 
number state. Note that this impulse function has a periodicity of 27r/u; due to the argument 
in the exponential functions and the 7 t/lu periodicity of the squared sine argument in Tricomi 
confluent hypergeometric function. This is shown in Fig. [6] for an initial field impinging 
the m — 24 waveguide in a Heisenberg-Weyl lattice with parameter A = 0.4a;. A displaced 
number state is nothing else than the state obtained by the action of Glauber displacement 
operator over a Fock state, which may be the origin of the Glauber-Fock lattice moniker 
given to specific realizations of lattices with underlying Heisenberg-Weyl symmetry [231126] . 
Note that these well known results for lu = 0 are recovered by taking the limit of the impulse 
function, 


ISHOA*) = z), (51) 

a;—>0 ’ 

= — = e~^ A ^ 2 ( i\z) n ~ m U (—m, n — m + 1, |A^| 2 ) , (52) 

ym\n\ 

and in this particular realization the periodicity is lost as shown in Fig. [7] A more com¬ 
plicated realization belonging to this symmetry class is provided by what we will call a 
Lewis-Ermakov lattice, where first and second neighbors couple and effective refractive in¬ 
dices and couplings may include propagation distance dependence m- 


V. CONCLUSION 

We have used the well known Gilmore-Perelomov coherent states from quantum me¬ 
chanics [T6HT8] to provide closed form analytic dispersion relations, normal modes and im¬ 
pulse functions for three types of photonic lattices with underlying S77(2), S77(l, 1) and 
Heisenberg-Weyl symmetries. This symmetry based approach allowed us to analyzed the 
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n 

100 


FIG. 7. Intensity at the n-th waveguide, A, z )| 2 , for light initially impinging the m = 24 

waveguide of a Heisenberg-Weyl array with parameter ui = 0. 

characteristics of these photonic lattices straight from the impulse functions; in particular, 
their periodicity or lack of it. We want to emphasize the fact that the propagator for any 
given tight-binding array of photonic waveguides may be calculated in this way as long as 
the dynamical group is found, even in the case of arrays were the parameters depend on the 
propagation distance HUES]. Furthermore, these photonic lattices can be seen as optical 
simulators of a class of Gilmore-Perelomov generalized coherent states. 
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